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Results of classical force field geometry optimizations for twisted graphene nanoribbons with 
^ a number of twists Nt varying from to 7 (the case Nt=l corresponds to a half-twist Mobius 

' O nanoribbon) are presented in this work. Their structural stability was investigated using the Brenner 

O reactive force field. The best classical molecular geometries were used as input for semiempirical 

calculations, from which the electronic properties (energy levels, HOMO, LUMO orbitals) were 
^ computed for each structure. CI wavefunctions were also calculated in the complete active space 

o 

OO framework taking into account eigenstates from HOMO-4 to LUMO+4, as well as the oscillator 

o 

strengths corresponding to the first optical transitions in the UV-VIS range. The lowest energy 
molecules were found less symmetric than initial configurations, and the HOMO-LUMO energy 
gaps are larger than the value found for the nanographene used to build them due to electronic 



localization effects created by the twisting. A high number of twists leads to a sharp increase of the 
HOMO — > LUMO transition energy. We suggest that some twisted nanoribbons could form crystals 
stabilized by dipolar interactions. 
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I. INTRODUCTION 

In 1858, the German mathematician August Ferdinand Mobius investigated the properties of a surface obtained 
by gluing the extremities of a twisted band. More than one century later, molecules with the shape of a half-twist 
Mobius band were synthesized for the first time [1 . In 1964, Heilbronner [2 proposed the definition of a Mobius 
aromaticity for cyclic molecules, foreseeing that Mobius aromatic hydrocarbons with 4n p-electrons should be stable, 
in contrast with the Hiickel rule prediction on the instability of antiaromatic hyrdrocarbons with An p-electrons. 
According to Heilbronner calculations, the twist changes the electronic structure of a molecular ring to a closed shell 
configuration. It was also demonstrated that the topological change to a half-twist Mobius band leads to a simple and 
general relationship between the electronic energy eigenvalues of the corresponding molecules [3], giving an elegant 
description of pure p-like systems over a nonorientable (one-sided) surface. 

In a more recent work, Mauksch et al. |4] presented a computational reinterpretation of experimental data showing 
that (Cll)tj" could be a Mobius aromatic cyclic annulene with An tt electrons. Electronic properties of ring compounds 
were discussed theoretically in connection with Mobius aromatic properties by some authors O E] |7l [H [9] . Martfn- 
Santamaria and Rzepa j5j presented an analysis of the main features of Mobius annulenes considering the 7r-molecular 
orbital correlation between the planar Hiickel configuration and the twisted C2 symmetric Mobius system. Their 
calculations were performed within both the frameworks of the semiempirical SCF-MO method and the restricted 
B3LYP/6-31G(d) level of ah initio Density Functional Theory (DFT). They observed that the form of the calculated 
AMI and B3LYP/6-31G(d) orbitals was very similar, with minor variations in the relative order of the occupied 
orbitals. They concluded that when ring C-H bonds are present in a cyclic annulene the Mobius orbitals can interact 
significantly with the a framework, creating a "ghost" 7r-like Hiickel orbital. In other paper [6J, Martin-Santamaria 
and Rzepa obtained AMI and ah initio geometries of Mobius bands formed by imposing one, two or three twists to 
cyclacenes of various lengths, showing that there is a localisation of the twist over 2-4 benzene rings, and reporting 
the geometries, highest occupied molecular orbitals, electrostatic potential and charge distributions for neutral, 6-1- 
and 6- forms of these molecules. Rzepa [7] carried out B3LYP and KMLYP/6-31G(d) calculations that predicted a 
double-helical and chiral conformation of [14]annulene with the topological properties of a double-twist Mobius band. 
In another work [9 , the same author summarized the evidence to support Heilbronner's prediction of a counterpart 
to planar Hiickel aromatic rings with Mobius topology and proposes higher order twisted systems (paradromic rings 
or Listing rings) as interesting candidates for future study. 

The synthesis of a stable Mobius aromatic hydrocarbon was obtained for the first time in 2003 by Ajami et al. [TU] 
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who combined a 'normal' aromatic structure (like benzene, with trigonal planar sp^-hybridized atoms) and a 'belt- 
like' aromatic structure (like the surface of a carbon nanotube, with pyramidalized sp^ atoms), creating a Mobius 
compound stabilized by an extended 7r-system. X-ray characterization of crystals grown from the synthesized Mobius 
hydrocarbons showed good agreement with DFT calculations. C2 symmetric Mobius compounds lead to red crystals, 
whereas the Cs symmetric Hiickel counterpart compound is colourless. Observed trends in bond-length equalization 
and stabilization energy point to a moderately aromatic character of the Mobius structure, whereas the corresponding 
Hiickel structure is non-aromatic. Further Mobius isomers were presented by Ajami et al. in 2006 [llj with details of 
their preparation as stable compounds and an investigation on their energy, geometry and magnetic parameters. We 
note also that Starostin and Van Der Heijden |12j have used the invariant variational bicomplex formalism to find the 
first equilibrium equations for a wide strip subjected to large deformations, which can be useful to understand the 
relationship between geometry and physical properties of nanometric Mobius structures. A review of the most recent 
developments on the design of Mobius molecules can be found in reference jlB] . 

The aim of this work is to investigate theoretically the structural stability, optimal geometry, and the electronic and 
optical properties of a class of twisted bands built from graphene nanoribbons, including the case of a half-twist Mobius 
graphene strip. Indeed, graphene is a novel material with very interesting and promising properties, such as thermal 
conductivity, hardness, and extraordinary electronic properties, like being a zero gap semiconductor with a Dirac-like 
energy spectrum close to the Fermi level (or, in other words, electrons in this region have zero effective mass and are 
analogous to photons in Einstein's theory of relativity) [TH [TSl HSJ EZl UHl EHl UHl 1211 122 • Zigzag nanoribbons from 
graphene have been investigated using first principles calculations and it was shown that they can carry a spin current 
[23j . Under Mobius boundary conditions, these nanoribbons have a magnetic domain wall whose width depends on 
the Coulomb interaction [53] . A first principles study on the spin states of zigzag graphene nanoribbons has found 
a triplet ground state for a Mobius cyclacene and an open-shell singlet ground state for the corresponding two-edge 
counterpart. Wider nanoribbons with Mobius topology exhibit an increasing total magnetization with the ribbon 
length 

II. COMPUTATIONAL DETAILS 

A single rectangular nanographene (see Figure 1(a)) with armchair and zigzag sides was used for all calculations 
we have carried out. We can describe it by using two integer parameters: its length (L) and width {W). The length 
L is obtained by counting the number of C-C dimers along (and parallel to) the armchair edge. The width W , on the 
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other hand, is given by counting the number of C-C dimers orthogonal to the zigzag border. Thus our nanographene 
has i = 29 6.8 nm) and = 5 0.5 nm), its molecular formula being CisoHgo (we have chosen the values for 
L and W in such a way to achieve a round number of carbon atoms). Closed nanoribbons can be built by joining 
the nanographene zigzag ends. If we twist the nanographene along its length by an angle of 180° and join the zigzag 
extremities, we obtain a half-twist Mobius nanoribbon (Figure 1(c)). If we twist the nanographene by an angle of 
N( X 180" and join the zigzag extremities, we build a A^t-twisted nanoribbon (Figure l(d)-(i)). In this work, we have 
chosen structures with Nt varying from to 7. Dangling bonds were passivated with hydrogen atoms during the 
simulations. 

We started out by performing classical geometry optimizations for each nanoribbon using the Universal force 
field, available in the Forcite code. The following convergence tolerances were adopted: 2 x 10^^ kcal/mol (total 
energy variation), 0.01 kcal/mol/nm (maximum force per atom), and 10^^ nm (maximum atomic displacement). An 
algorithm using a cascade of the steepest descent, adjusted basis set Newton-Raphson and quasi-Newton methods 
was selected. Non-bond energies (van der Waals and electrostatic) were calculated using an atom based summation 
method with a cutoff distance of 6 nm. After this optimization, an annealing simulation was carried out to search the 
conformational space for lowest energy structures through cycles of increasing-decreasing temperatures of a classical 
dynamics in order to prevent conformations that correspond to local energy minima with higher energies. For the 
annealing, we have use the following parameters: NVE ensemble, time step of 1 fs, 100 annealing cycles, initial (final) 
temperature of 300 K (800 K), 50 heating ramps per cycle and 100 dynamic steps per ramp. After each anneal cycle, 
the lowest energy structure was optimized following the same criteria of the first geometry optimization. Aiming to 
evaluate the structural stability of the twisted nanoribbons, a classical molecular dynamics using the Brenner reactive 
forcefield [571 HH] and the GULP code was also performed starting from the optimized geometries and using a NVE 
ensemble with a temperature of 500 K, time step of 1 fs, equilibration time of 500 ps, and production time of 500 ps. 

The electronic structure was simulated for all nanoribbons using the AMI semiempirical Hamiltonian. Previous 
reports in the literature [SOJ EH suggest that both AMI and first principles Hartree-Fock results are really 
comparable for both carbon and silicon cages and provide firm support to the applicability of the AMI method to study 
all-silicon and all-carbon clusters. Thus we think it is appropriate to use the semiempirical AMI Hamiltonian to study 
the graphene nanoribbons as well. Wavefunctions were calculated within the restricted Hartree-Fock approximation 
(ground state in a singlet spin state), and a geometry optimization was carried out starting from the best molecular 
conformations achieved through the classical annealing. The lowest vibrational frequencies were monitored for negative 
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values to ensure that the semiempirical conformations correspond to local energy minima. Electronic energies and 660 
self-consistent molecular orbitals (330 occupied) were then obtained. After the semiempirical geometry optimization, 
the configuration interaction (CI) wavefunctions for the singlet ground state and the first 63 excited states were 
calculated using the complete active space method and taking into account the molecular orbitals (MOs) from HOMO- 
4 to LUMO+4 (10 MOs). From these CI wavefunctions, the energies corresponding to the first optical transitions 
and their respective oscillator strengths were computed. 

III. RESULTS AND DISCUSSION 
A. Geometry and stability of the twisted nanoribbons 

In Table 1 we present the space point groups of symmetry for the nanoribbons and the rectangular nanographene 
used to build them. The nanographene has I?2/i symmetry (8 symmetry operations), whereas the initial configuration 
for the 0-nanoribbon is the most symmetrical of all, belonging to the Di^h group (60 symmetry operations). The 
classically optimized 0-nanoribbon (after the annealing) exhibits a much smaller degree of symmetry (C2 point group). 
The semiempirical optimization, on the other hand, leads to a 0-nanoribbon with symmetry belonging to the 
point group (10 symmetry operations). Much of this symmetry breaking is due to the insertion of hydrogen atoms 
to passivate the dangling bonds at the nanoribbon edges, which creates more degrees of freedom for the annealing 
simulation to explore asymmetrical geometries. Indeed, the symmetry decreases with respect to the initial inputs for 
all the other nanoribbons, from Nt = I (half-twist Mobius) to 7, after the annealing and semiempirical simulations. We 
can say therefore that the twist localization created by such symmetry breaking is not merely a quantum mechanical 
efi^ect, but a more general mechanical feature of twisted ribbons, as pointed in the paper of Starostin and Van Der 
Heijden JTT. The half-twist Mobius nanoribbon {Nt = 1) has an initial C2 point group, but ends at a Ci geometry 
for both the classical and semiempirical results. Nanoribbons with Nt = 2,4,7 start with C2 symmetry and end 
with Ci geometries as well. The Nt — 3,6 molecules end both at C2 symmetry in the classical and Ci in the 
semiempirical simulation. Finally, for Nt — 5, we have an initial atomic structure with C5 symmetry that ends with 
C2 point group for both classical and semiempirical calculations. In comparison, Martfn-Santamarfa and Szepa [S], 
using semiempirical and ab initio methods, have studied single, double and triple twisted Mobius cyclacenes, and 
have shown that the optimized geometries have localized twists and reduced symmetry, agreeing with the results we 
obtained for the twisted armchair nanoribbons. We believe this agreement points to localizing behavior as a general 
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topological feature of Mobius systems. 

In a classical molecular mechanics simulation the total energy of a system is usually decomposed in contributions 
due to the stretching of bonds, the bending of angles between bonds from the same atom, the torsion of dihedral 
angles, the existence of out-of-plane interactions, and the non-bond interactions (van der Waals, Coulomb, hydrogen 
bonds). Table 2 shows the total energy and the contribution from valence and non-bond terms (only van der Waals) 
for all nanoribbons from Nt — to 7. We can see that the total energy increases with the number of twists Nf. A 
quadratic fit shows that the total energy i?TOT as a function of Nt (with maximum error of about 5%) can be given 
by ETOT{Nt) = 311.94 - 16.82iVf -I- 11A7N^. The Nt = 7 strip has a total energy 2.6 times larger than the total 
energy for a nanoribbon with Nt = 0. The variation of energy from Nt = to Nt — I (half-twist Mobius) is 24.91 
kcal/mol, corresponding to a relative increase of about 8.4%. This total energy increase behavior was also observed 
by Fowler and Rzepa using a pure Hiickel treatment for cycles with an arbitrary number of half-twists [33j . 

Looking at the terms contributing to the total energy, we see that the energy related to the stretching of bonds is 
practically the same (around 30 kcal/mol) for all nanoribbons from A't = to 5, indicating that the change in bond 
lengths is not remarkable for structures with Nt < 6. On the other hand, for Nt — 6,7 we note a more significant 
increase of the bond energy term, pointing to an increase in the bond strain as the number of torsions crosses a 
structural threshold. Bond energies contribute with 11% of the total energy at most (for Nt — 0) and this figure is 
smaller for larger values of Nt (reaching about 5% for Nt ^ 7). The energy term related to the bending of covalent 
bonds, however, is increasingly more important as we switch from to 7-nanoribbons. It is smaller than the bond 
energy for Nt < 4, but increases to approximately 12.5% of the total energy for Nt ~ 7. 

It is expected that, as the nanoribbons are subjected to more and more twists, the energy related to the torsion of 
dihedral angles must increase almost linearly with Nt- A linear interpolation, however, is not a good approximation 
for the behavior observed in our calculations. As occurred with the total energy, the torsion energy increase with 
Nt can be modeled by a parabolic fit according to i?TORSiON(A''f ) = 65.9 -I- 11.3iVf -I- 6.0A^(^ (maximum error of 9%). 
The relative contribution of the torsion term to the total energy increases from 21.3% {Nt — 0) to 57.3% {Nt = 7), 
becoming dominant (larger than the van der Waals energy) for Nt > 2. The van der Waals is the only non-bond 
energy taken into account in the classical calculations carried out in this work, and varies from 61.6% {Nt ~ 0) to 
approximately 18.5% {Nt — 7) relative to the total energy. Its absolute value decreases from « 183 kcal/mol {Nt — 0) 
to K, 130 kcal/mol {Nt — 5), and the increases to w 137 kcal/mol {Nt — 6), and k, 142 kcal/mol {Nt — 7). Finally, 
the inversion term is the smallest from all contributions for iVj < 6, ranging from only 1.79 kcal/mol {Nt — 0) to 
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53.18 kcal/mol {Nt = 7). Notwithstanding, it is the contribution from the inversion energy to the total energy that 
presents the largest relative variation of all terms considered in the calculations, varying from 0.6 % {N^ — 0) to 6.9 
% iNt=7). 

The chemical stability of all twisted nanoribbons was also investigated through a classical molecular dynamics 
simulation using the Brenner forcefield, which provides a reactive bond-order potential useful to model organic systems 
[271 12H]- The simulation temperature was set to 500 K in a NVE ensemble, with a time step of 1 fs and a total time 
of 1 ns. We have not observed any bond dissociation for all nanoribbons studied in this work, which indicates they 
are chemically stable under the simulated conditions. Beyond that, we mention also the recent work by Rzepa |34j 
showing that the stabilization of twisted moieties could be associated with the reduction of local twists of adjacent 
p^-p^ orbitals by the coiling of the central axis of the ring into three dimensions through writhing. Such coiling can 
indeed be observed in the optimized structures of this work, as shown in Figs. 1 (results from classical simulations) 
and 3 (quantum semiempirical results) . 

B. Electronic structure 

As observed for the classical molecular mechanics simulations, the total energy for the geometries optimized using 
the semiempirical approach shows a quadratic dependence on the number of twists. The heat of formation {HF) of 
the nanoribbons can be given with good accuracy (relative error always smaller than 3%, being maximal for Nf = 1, 
the Mobius band) by HF{Nt) = 784.58 + 3.43A'^t + 15MN^ kcal/mol. One sees in Figure 2 the highest occupied 
molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) eigenenergies for all nanoribbons (and 
the classically optimized nanographene) , as well as the HOMO-LUMO energy gap, the CI energy gap and the energy 
of the first allowed optical transition. The HOMO and LUMO eigenenergies behave in opposite ways as Nt increase, 
with the HOMO (LUMO) tending to increase (decrease) slightly with some oscillation as Nt varies from to 5 and 
decreasing (increasing) sharply for Nt — 6,7. The HOMO-LUMO gap Eg, defined as the difference between the 
HOMO and LUMO energies, Eg — E'lumo ^ Eb_omo, behaves qualitatively like the LUMO as a function of Nt, 
starting from 3.67 eV for the nanographene rectangle, increasing to 4.65 eV for the 0-nanoribbon and decreasing to 
4.58 eV for the 1-nanoribbon (with half-twist Mobius topology). A comparison with previous calculations within 
the framework of density functional theory in the generalized gradient approximation (DFT-GGA, PBE exchange- 
correlation functional) for finite armchair nanographenes [35^ shows that, in the case of a (W = 5, L = 12) molecule, 
the HOMO-LUMO energy gap is 0.6 eV and, for an infinite 1D-(W^ = 5) nanographene Eg,K, 0.3 eV. Considering 
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that the graphene rectangle in our work has L = 29, we could expect a DFT-GGA energy gap possibly between 0.3 
eV and 0.6 eV. To check for this, DFT simulations were carried out (using only the classically optimized geometry 
of the flat nanographene) and we obtained the same value of « 0.11 eV for the HOMO-LUMO energy gap using the 
GGA-PBE, GGA-BLYP, and LDA-VWN exchange-correlation functionals, values smaller than we could expect based 
in the prediction for the ID case of reference [35.,. It seems that the HOMO-LUMO gap of the graphene finite ribbon 
with = 5 as a function of L does not decrease monotonically, but oscillates about the value for the L — oo case. 
Comparing the DFT estimates with the semiempirical result (3.67 eV), we see that it predicts a much higher HOMO- 
LUMO gap. We must remember here that, despite including some electron correlation energy via its parameters, the 
AMI Hamiltonian produces energy gaps close to the ones obtained through the Hartree-Fock approximation [301136] . 
which does not take into account any correlation effects and, therefore, widens the energy gaps in comparison to 
experimental values. DFT methods, on the other hand, tend to overestimate electronic correlation, underestimating 
energy gaps. So we must exercise some caution when analyzing the semiempirical and DFT results. As we will 
see, the CI calculations include some improvements on estimating electronic correlation energy, predicting energy 
gaps between the ground state and the first excited state intermediary between the DFT and AMI estimations (and 
hopefully more accurate than both) for all nanoribbons. For Nt ~ 2, the AMI HOMO-LUMO energy gap increases 
in comparison with the Nt = \ (half-twist Mobius) case to 4.66 eV and then decreases for Nt = i down to 4.56 eV. 
The Nt = A gap is very close to Eg{Nt = 3), being 4.55 eV. The smallest gap is obtained for the iVf = 5 molecule, 
Eg{Nt = 5) = 4.42 eV. Finally, Eg{Nt = 6) = 5.20 eV and Eg{Nt = 7) = 5.45 eV (the largest ones). 

To explain why the HOMO-LUMO energy gap increases sharply for Nt > 5, we have plotted the HOMO and 
LUMO orbitals for all nanoribbons and the original nanographene rectangle, as shown in Figure 3. One can see 
that the nanographene has hydrogen atoms passivating the dangling bonds in its extremities which are absent in 
the nanoribbons. The nanographene HOMO (LUMO) orbital is delocalized over its length, resembling a metallic 
quantum state. Indeed, there is a rule for ID (infinite) nanographenes according to which these nanostructures are 
metallic ii W — 3N — 1, where N is an integer [371 |3H1 [351 HO]- As the finite (OD) nanographene studied here has 
M^ = 5 = 3x2— 1, one can see that in this particular case such rule holds. The phase of the HOMO orbital alternates 
its sign when we move along the nanographene length. We will denote this amplitude pattern 7ri. On the other hand, 
the LUMO orbital has an amplitude that changes its sign as we move across the width of the nanographene, thus we 
have chosen to indicate this pattern using the ttW abbreviation. Both patterns were also observed in the previously 
mentioned calculations using DFT-GGA and DFT-LDA methods to check the nanographene HOMO-LUMO energy 
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gap. Nanoribbons with Nt varying from to 5 have HOMO (LUMO) ampUtude phase with wW {ttL) aspect, swapping 
the order of patterns observed for the original nanographene. For A^^ = 3, 4, 5, electrons occupying the HOMO and 
LUMO states clearly present a pronounced degree of localization (Figure 3(e,f,g,e',f',g')). For Nt = 6,7, one can see 
amplitude phase structures with the same order observed for the nanographene, and for = 7 electrons in HOMO or 
LUMO states become strongly localized. So we conclude that the formation of a closed nanoribbon initially inverts the 
phase patterns of frontier orbitals observed for the armchair nanographene, but the twisting localizes the electrons and, 
ultimately, restores the order of the nanographene phase patterns. The electronic localization contributes to increase 
the HOMO-LUMO energy gap. When iVj = 6, 7, the restoring of the phase-patterns produces the sharp change noted 
in all plots displayed in Figure 2. At last, besides the HOMO and LUMO orbitals, we have also investigated the lowest 
energy 7r-like orbital for each Nt and observed that, for even values of Nt, there is a continuous orbital ribbon on each 
face of the strip (even Nt implies in an orientable surface) with no phase shifts, whereas for odd Nt the corresponding 
orbital ribbons have at least one phase shift, as shown for A^^ = 1, 2 in Figure 4. 

According to the CI calculations, the singlet ground state and the first excited state - which is a triplet for all 
nanoribbons and the nanographene - are separated by an energy gap Eqi about 3 eV smaller (on average) than the 
HOMO-LUMO gap calculated using the AMI Hamiltonian. Such energy gap reduction is due to the incorporation 
of some electronic correlation energy in the CI method. As Nt is switched from to 7, Eci{Nt) varies qualitatively 
much like Eg. For the nanographene, Eqi = 1.16 eV. The 0-nanoribbon has Eci = 1.86 eV and the half- twist Mobius 
band has Eci = 1.75 eV. For the nanoribbons up to Nt = 6 we have Eci{Nt = 2) = 1.81 eV, Eci{Nt = 3) = 1.59 
eV, Eci{Nt = 4) = 1.63 eV, Eci{Nt = 5) = 1.38 eV (the smallest of all energies), and Eci{Nt = 6) = 1.93 eV. It 
can be seen that Eci{Nt) oscillates: an odd Nt nanoribbon has a smaller energy gap than the value predicted for the 
even Nt — 1 one. Exception to this rule is Eci{Nt = 7) = 2.03 eV. Optical transitions between singlet and triplet 
states are forbidden by the spin selection rule, so the first excited state accessible through the absorption of photons 
has an energy larger than the first triplet excited quantum level. The plot at the bottom right of Figure 2 shows the 
first optical absorption transition energy ascribed to the nanoribbons, EopTi^t)- It follows closely the plot for Eci- 
The absorbed wavelengths corresponding to such transitions correspond to infrared (nanographene, Nt = 1,3,4,5), 
orange {Nt = 0), red {Nt = 2), yellow {Nt = 6), and green light {Nt = 7). The intensity of the absorption peaks is 
related to the oscillator strength /, shown at the bottom of Figure 4. One can see that the nanographene has a very 
strong oscillator strength in comparison with the closed nanoribbons. f{Nt) oscillates as we move from A^^ = to 
Nt — 7. For Nt = 0,2,4,6 (even Nt) we have the minima of f{Nt) whereas A^^ = 1,3,5,7 (odd Nt) correspond to 
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maxima. The Nt = ribbon has the smallest value of /, w 10~^ a.u., due to the symmetries of the ground state and 
first singlet excited state (indeed, one can say that this transition is practically forbidden). The first optical transition 
with significant oscillator strength for the 0-nanoribbon involves a singlet state with energy of 2.53 eV (corresponding 
to blue light), with / w 9.4 a.u. For Nt — 1 (half-twist Mobius ribbon), the largest value of / occurs for Ed = 2.43 
eV (blue light), and for Nt — 6,7, f is maximum for Eci = 2.86 eV and Eqi — 2.93 eV (both equivalent to a blue 
wavelength), respectively. f{Nt — 5) (« 0.8 a.u.) is the largest of all first optically active electronic transitions shown 
in Figure 2. 

Finally, at the top of Figure 5 the magnitude of the ground state dipole moment p as a function of Nt is shown. 
It was demonstrated that the use of configuration interaction including only HOMO and LUMO in the space of 
configurations leads to AMl+CI values for p with practically the same quality as first principles results 41J. The 
Nt — 3,6 nanoribbons have the largest (and very close) dipole moments, « 0.74 D. In contrast, the Nt ~ ribbon 
has the smallest value, ~ 10^^ D (5 times smaller than the dipole moment of the nanographene rectangle). In the 
case of a half-twist Mobius strip, p{Nt) = 0.32 D. p{Nt) has an order of alternate maxima and minima for < A^t < 5 
(minima for and even Nt, maxima for odd Nt). This order is broken down for Nt — 6. We think that crystals 
made from some of the nanoribbons (specifically those with Nt odd and the Nt = 6 case) could be stabilized through 
dipole-dipole interactions, as it occurs with many organic compounds. 

IV. CONCLUSIONS 

In this paper, we presented the results of computational simulations to obtain structural and electronic properties of 
twisted nanoribbons built from a nanographene strip. Optimized geometries were found through classical molecular 
mechanics calculations in three steps: (i) first geometry optimization starting from symmetric initial structures; 
(ii) thermal annealing to search for geometries with lower total energy; (iii) second geometry optimization of the 
nanostructures with lowest total energy found after the annealing. A stability study was performed on the optimal 
geometries to check for the dissociation of chemical bonds. All nanoribbons were stable for a temperature of 500 K 
and simulation time of 1 ns. The total energy as a function of the number of twists was fitted to a parabolic curve 
with good accuracy, the largest fit error being observed for the Nt — I (half-twist Mobius band) case. Optimized 
geometries have lower symmetry in comparison with starting configurations. 

The electronic structure was calculated for all nanoribbons using the AMI semiempirical Hamiltonian and a CI 
calculation was carried out to improve the estimate of the HOMO-LUMO energy gap, which is larger for the twisted 
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molecules in comparison to the flat metallic nanographene. We noted a sharp increase of the HOMO— >LUMO 
transition energy as the number of twists varies from Nf = 5 to = 6,7, indicating that there is a change of 
electronic structure related to the increase of torsion in the twisted nanoribbons. When a closed nanoribbon {Nt = 0) 
is formed, the HOMO-LUMO wavefunction appearance is inverted in comparison to the armchair nanographene case. 
Increasing Nt leads to electronic localization in the frontier orbitals and restores the nanoribbon HOMO-LUMO 
patterns to those of the nanographene. AMl+CI energy gaps are much smaller than the ones calculated using the 
AMI Hamiltonian, but probably much larger than the DFT-GGA estimates. The nanographene oscillator strength / 
for the first optically active electronic transition is stronger than the corresponding ones for the closed nanoribbons. 
For Nt = 0,2,4,6 (even Nt), f{Nt) has minima, whereas Nt = 1,3,5,7 (odd Nt) correspond to maxima of /. The 
first allowed optical transition for the case Nt = between the singlet ground state and the singlet excited state 
is forbidden due to the wavefunction symmetries. For Nt = I (half-twist Mobius nanoribbon), / is maximum for 
Eci = 2.43 eV. Crystals made from the nanoribbons with Nt odd and Nt = 6 could possibly be stabilized through 
dipolar interactions. 
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TABLE I: Space point groups of the graphenc nanoribbons studied in this work: initial geometry (second column), classically 
optimized geometry (after annealing, third column), and semiempirical optimized geometry (fourth column). The second row 
shows the space point group of the nanographene rectangle used to build the nanoribbons. 



Number of 


Initial 


Optimized 


Optimized 


twists {Nt) 


geometry (annealing) (semiempirical) 


Nanographene 


D2h 









Dl5h 


Ca 


Csji 


1 


C2 


Ci 


Ci 


2 


C2 


Ci 


Ci 


3 


Cs 


C2 


Ci 


4 


C2 


Ci 


Ci 


5 






C2 


6 


C3 


C2 


Ci 


7 


C2 


Ci 


Ci 



15 



TABLE II: Total energy and its components (bond, angle, torsion, inversion, and van der Waals) for the optimized molecules 
after the classical annealing. The percentual of total energy for each contribution is also shown. All energies are given in 
kcal/mol. 





Total 


Bond 


Angle 


Torsion 


Inversion 


vdW 


Nt 


Energy 


Energy 


Energy 


Energy 


Energy 


Energy 





296.74 


31.97 


16.99 


63.22 


1.79 


182.77 


1 


321.65 


30.18 


18.31 


91.78 


3.57 


177.81 


2 


331.70 


29.55 


18.00 


103.58 


4.47 


176.10 


3 


372.52 


28.97 


22.14 


158.38 


10.67 


152.36 


4 


423.25 


29.32 


30.42 


207.13 


15.90 


140.48 


5 


501.19 


30.74 


42.80 


273.71 


24.32 


129.62 


6 


613.12 


33.58 


54.99 


351.98 


35.50 


137.07 


7 


770.37 


37.66 


96.03 


441.31 


53.18 


142.19 






Bond/ 


Angle/ 


Torsion/ 


Inversion/ 


vdW/ 


Nt 




Total (%) Total (%) Total (%) Total (%) Total (%) 







10.77 


5.73 


21.30 


0.60 


61.60 


1 




9.38 


5.69 


28.53 


1.11 


55.29 


2 




8.91 


5.43 


31.23 


1.35 


53.08 


3 




7.78 


5.94 


42.52 


2.86 


40.90 


4 




6.93 


7.19 


48.94 


3.76 


33.18 


5 




6.13 


8.54 


54.61 


4.85 


25.87 


6 




5.48 


8.97 


57.41 


5.79 


22.35 


7 




4.89 


12.47 


57.29 


6.90 


18.45 




FIG. 1: (a) Rectangular nanographene used to build the twisted nanoribbons. The twisted nanoribbons, after classical annealing 
and geometry optimization, are presented; (b) Nt = nanoribbon; (c) Nt = 1 nanoribbon (half- twist Mobius strip); (d) Nt = 2; 
(e) Nt = 3; (f) = 4; (g) = 5; (h) Nt = 6; (i) N* = 7. 
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FIG. 2: HOMO and LUMO energy levels, HOMO-LUMO energy gap, CI energy gap (-Eci, solid squares) and the energy of the 
first allowed optical transitions (-Eqpt, open squares) for all the nanoribbons and the nanographene rectangle (NG). 




FIG. 3: HOMO and LUMO orbitals for the: (a,a') Rectangular nanographene; (b,b') Nt = nanoribbon; (c,c') Nt = 1 
nanoribbon (half-twist Mobius strip); (d,d') Nt = 2; (e,e') Nt = 3; (f.f) Nt = 4; (g,g') Nt = 5; (h,h') Nt = 6; (i,i') Nt = 7. 
Isosurfaces correspond to a wavefunction amplitude of 0.002. 




FIG. 4: Lowest energy 7r-like orbitals for Nt — 1 (a) and Nt = 2 nanoribbons. For the Nt 
twist the appearance of a phase shift in the wavefunction. There is no phase shift in the Nt 



— 1 case, we note at the locahzed 
= 2 ribbon. 
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FIG. 5: Dipole moment and oscillator strength of the first allowed optical transition for each graphene nanoribbon and for the 
nanographene molecule (NG). 



